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We present a general scheme for performing a simulation 
of the dynamics of one quantum system using another. This 
scheme is used to experimentally simulate the dynamics of 
truncated quantum harmonic and anharmonic oscillators us- 
ing nuclear magnetic resonance. We believe this to be the 
first explicit physical realization of such a simulation. 

PACS numbers: 03.67.-a,76.60.-k 



In 1982, Richard Feynman proposed that a quantum 
system would be more efficiently simulated by a computer 
based on the principles of quantum mechanics rather 
than by one based on classical mechanics Recently, 
it has been pointed out that it should be possible to effi- 
ciently approximate any desired Hamiltonian within the 
standard model of a quantum computer by a sparsely 
coupled array of two-state systems Many of the 

concepts of quantum simulation are implicit in the av- 
erage Hamiltonian theory developed by Waugh and col- 
leagues to design NMR pulse sequences which implement 
a specific desired effective NMR Hamiltonian [|] . Here we 
show the first explicit simulation of one quantum system 
by another; namely the simulation of the kinematics and 
dynamics of a truncated quantum oscillator by an NMR 
quantum information processor |^,^. Quantum simula- 
tions are shown for both an undrivcn harmonic oscillator 
and a driven anharmonic oscillator. 

A general scheme for quantum simulation is summa- 
rized by the following diagram: 
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The object is to simulate the effect of the evolution 

\s) \s{T)) using the physical system P. To do this, 
S is related to P by an invertible map (j) which deter- 
mines a correspondence between all the operators and 
states of S and of P. In particular, the propagator U 
maps to Vt = (j)~^U(j). The challenge is to implement Vt 
using propagators Vi arising from the available external 
interactions with intervening periods of natural evolu- 



sufficient class of simple operations (logic gates) are im- 
plementable in the physical system, the Universal Com- 
putation Theorem ^-p^ guarantees that any operator 
(in particular Vt) can be composed of natural evolutions 
in P and external interactions. For unitary maps </>, we 
may write Vt = e~''^p^/'' where Tip = (j)~^Tis4> can be 
identified with the average Hamiltonian of Waugh. Af- 
ter \p) --^ \pt)i the final map (j)^^ takes \pt) l^iT)) 
thereby effecting the simulation \s) \s{T)). Note that 
Tig (T) can be a time dependent Hamiltonian and that T 
is viewed as a parameter when mapped to P. This im- 
plies that the physical times ti{T) are parameterized by 
the simulated time T. 

Liquid state NMR quantum computers are well suited 
for quantum simulations because they have long spin re- 
laxation times (Ti and T2) as well as the flexibility of 
using a variety of molecular samples. In particular, the 
coupling between the nuclear spins, usually dominated 
by the 'scalar' coupling (J), may be reduced at will by 
means of radiofrequency pulses. Typically spin 1/2 nuclei 
are used. Thus, the kinematics of any 2^ level quantum 
system could be simulated using a given iV-spin molecule. 

We chose to simulate a quantum harmonic oscilla- 
tor (QHO) with a 4 level, 2-spin system P being the 
two proton nuclear spins in 2,3- dibromothiophene. The 
Hamiltonian of a QHO is Hqho = fin{N + ^) = 
^n^^{n + ^)\n){n\, where i7 is the oscillator frequency, 
and 1 71) are the orthonormal eigenstates of the number 
operator N. Since the nuclear spin eigenstate space is fi- 
nite dimensional (4 levels, in this case), only a truncated 
version of the infinite dimensional oscillator was simu- 
lated. However, as noted above, for N spins, there are 
2^ levels. A convenient unitary mapping, (j>, between the 
energy eigenstates of the QHO and a 2-spin system is: 
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While any of a number of mappings would suffice, this 
mapping is convenient since An = ±1 correspond to al- 
lowed transitions in P. This mapping generalizes to a 
Gray code. Also note that the spin basis, permuted un- 
der (f), is now not in order of increasing energy in P. 

When truncated, e-*'^QffoT/;i jg mapped onto the 2- 
spin system as follows: 
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['^\y ' " - [^I^J] - [-57r/6]|^+^ - [G] , where the magnetic 
field gradient [G] destroys ofF-diagonal terms in the den- 
sity matrix. The sequence (||) for Vt then leads to 
\pt){pt\- Since the simulated system should evolve co- 
herently according to the diff'erence in energy levels of 
the various superpositions, \pt){'Pt\ above expresses a 
2iir dependence. NMR experiments are sensitive only 
to transverse dipolar magnetization, corresponding to the 
boxed components in the density matrices above. Thus 
a final read pulse is needed to rotate the e^*^^^ elements 
into view. The result manifests itself as a 2Q,T oscillation 
of the spectral peak heights as a function of the indirect 
dimension T . 

The dynamics of the truncated QHO states |0), (|0) -I- 
i|2)), and (|0) + |1> + |2) + |3)) were simulated. Eigen- 
states like |0) do not evolve, as the simulation in Figure 
1(a) shows. Fig 1(b) shows the 2^T oscillations discussed 
above for |0) In both Figs 1(a) and 1(b), [7r/2]^ 

read pulses were used. For \s) = |0) + |1) + |2) + |3), 
mixtures of SIT and 3S1T oscillations can be observed in 
the spectra. For example, the operator |0)(1| corresponds 
to I t T) (T i I which is a transition of spin 2. Thus the 
amplitude of the spin 2 peak will oscillate at QT. In 
Fig 1(c) SIT peak oscillations (on spin 2) are recorded 
while Fig 1(d) shows a superposition of SIT and BSIT (on 
spin 1). Since the two-spin system (P) has no natural 
triple quantum coherences, the latter coherence is en- 
tirely simulated. For Figs 1(c) and 1(d) read pulses were 
not required. 

In general, scaling the above to include more levels will 
depend on the various couplings between the added spins. 
For larger spin systems certain couplings are small and 
therefore severely limit the time scale of the experiment. 
For the truncated QHO however, an effective hamiltonian 
that is free of all couplings results from mapping the en- 
ergy eigenstate |fc) to the spin eigenstate corresponding 
to the binary representation of k in contrast to the Gray 
coding: 



expH(^|0)(0| + i|l)(l|) 
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Using the Pauli matrices {ca,, Uy, a^} as a basis for the 
2-spin density matrices |ll| , we may write 



Vt 



'iHpT/n 



e^^[i{al{l + \al)-2)nT]. (2) 



Implementing the operator (^) on the 2-spin system thus 
constitutes a simulation of the truncated QHO. This is 
easily done by making various refocussing adjustments to 
the physical 2-spin propagator e~'^p*'/^, obtained from 



the natural Hamiltonian TiJJ 



uj{))a'^ + -KJa^al), where a;i^2/27r are the resonance fre- 
quencies of spins 1 and 2, (w2 — wi)/27r = 226 Hz, loqI^-k 
is the spectrometer frequency (~ 400 MHz), and J is a 
scalar coupling strength (5.7 Hz). The following on res- 
onance (wo = wi) pulse sequence implements Vt for the 
simulated period rJT: 



1+2 



[T1/2 + T2]. 



(3) 



7r|;, ' " represents a tt angle radiofrequency 



The symbol [7r]i+2 
pulse, oriented along the y direction, on spins 1 and 2 
(corresponding to the Vi) and [r] represents a delay dur- 
ing which the 2-spin propagator e~^'^p^ acts. The time 
intervals are given by ri = S7T[l/(7r J) — 2/(^^2 ~ "^i)] and 

T2 = 2rir/ (w2 — wi). 

The experimental procedure is illustrated using |s) = 
|0) -t- i|2) as follows: 
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The initial state |p) = | t T) + *l i i) ^ is easily pre- 
pared from the (pseudopure Q) state | t ?)• This in turn 
is produced from the thermal equilibrium state of 2,3- 



dibromothiophene by the sequence [7r/4 



1-1-2 



[1/4J] 



'Hr, 
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This may be implemented by removing all scalar cou- 
plings and scaling all chemical shifts; for instance by 
methods analogous to "chemical shift concertina" se- 
quences introduced by Waugh |Q. 

The Hamiltonian for an anharmonic oscillator, 
T^AHO = fii^[{N -f ^) -t- Ai(^ + ^ f], where is the an- 
harmonicity parameter. The energy difference AEm be- 
tween the mth and (m -|- l)st energy level is A.Em = 
fiil[2fj,{m -f 1) -I- 1]. Radiation at the frequency AEm/h 
will drive a selective transition between these levels. The 
Hamiltonian for this selectively driven anharmonic oscil- 
lator is Haho + ^^^Rilm) {m+l\ + \m+l) (m\) , where V,fi 
is the Rabi frequency. Using the map (^, the |0) <-+ |1) 
driven truncated Hamiltonian in particular maps to Tip = 

inn [fial - 4(4^* + l)al{l + ^a])] + \nnn alii + a^). 
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This is implemented on 2,3-dibromothiophene by the fol- 
lowing pulse sequence 

Vt - [ri/2] - [vr],;, - [n/2] - [37r/4]^ - [r^] - [7r/4]^ . 

(4) 

For /i = —2/9, and ^b. = —2/91], the time intervals 
are determined by nT/2TT = 9Jr2/(2V2) = (9/2)(m - 
Tia;2/27r), where m is an integer. The receiver was set at 
ujq/2it = uji/2iT — J/2. Note that the map (|^) does not 
map the physical eigenstates to the simulated eigenstates 
(dressed states) of the driven oscillator emphasizing that 
knowledge of the eigenvalues/states of the simulated sys- 
tem is not assumed. Experimental results are shown in 
Figure 2. 

In these studies, we have only considered unitary evo- 
lution and have explored the quantum dynamics for sys- 
tems without dissipation. The decoherence [ p3[ intrinsic 
in our physical system (characterized by the longitudinal 
and transverse magnetization relaxation times, Ti and 
T2) limits the time of the experiment. This then limits 
the number of periods that can be simulated. Since the 
experimental (t) and simulated (T) time scales need not 
be identified with each other, this can be interpreted as a 
restriction either on or on T. In Fig 2 the visible decay 
due to T2 relaxation clearly shows this limitation. While 
decoherence can be controlled in principle by error cor- 
rection |p^-|l6||, it would be difficult to utilize this in the 
weakly polarized physical system used here. Moreover, 
thermal equilibrium will not necessarily map to another 
configuration that is thermal. Decoherence, itself, may 
be simulated through specific non-unitary evolution; in 
NMR for example by magnetic field gradients . 

The aspects available in the simulation: controlled 
kinematics and dynamics, a driving field, and decoher- 
ence suggest a very general tool with which to study other 
systems. 
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FIG. 1. NMR peak signals from 2,3- dibromothiophene 
demonstrating a quantum simulation of a truncated harmonic 
oscillator as implemented by Vt in (^). The various initial 
states express oscillations according to the energy differences 
between the eigenstates involved. The solid lines are fits to 
theoretical expectations, a, Evolution of the initial (pseudop- 
ure) state |0), showing no oscillation, b, Evolution of the 
initial state |0) + j|2), showing 2f2 oscillations, c, Evolution 
of the initial state |0) + |1) + |2) + |3) showing the fl and d 
oscillations. (Each trace in c and d is actually a combination 
of f2 and 3Q oscillations.) 
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FIG. 2. NMR peak signals from 2,3- dibromothiophene 
demonstrating a quantum simulation of a driven, truncated 
anharmonic oscillator as implemented by Vt in (^). When 
the (0,1) transition is selectively driven, the initial state |0) 
(o) undergoes Rabi (fin) oscillations to the |1) state, whereas 
the |2) state (*) does not evolve under the simulated Hamil- 
tonian. The exponential decay due to natural decoherence in 
P is clear. 



